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ABSTRACT 

We present a Sub millimeter Array study in the 1.3 mm waveband of the 
NGC 7538 IRS 1-3 massive star-forming region. The brightest core in the mm 
continuum map, MM1, harbors the IRS 1 young O star. The core has a gas 
temperature of about 245 K and shows spatially unresolved emission in complex 
organic molecules, all typical of a hot molecular core. Toward MM1, redshifted 
absorption is seen in molecular lines with different energies above the ground 
state. This absorption probes inward motion of the dense gas toward the central 
young O star, and the estimated mass accretion rate reaches 10 -3 M yr _1 . Mul- 
tiple outflows are seen in the CO and 13 CO maps. The gas mass of 50 M & and 
mass outflow rate of 2.5 x 10 -3 M yr _1 measured in CO line wings are domi- 
nated by the MM1 outflow, which is most likely driven by a fast wide-angle wind. 
Apart from MM1, we discover eight new dusty cores, MM2-9, within a projected 
distance of 0.35 pc from MM1. These cores show no counterpart in infrared or ra- 
dio continuum emission, while seven of them appear to be forming intermediate- 
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to high-mass stars. This manifests a deeply embedded star-forming component 
of the parent cloud of IRS 1-3. Apparently we are observing a Trapezium system 
in formation, and the system is presumably surrounded by a cluster of lower mass 
stars. 

Subject headings: ISM: individual (NGC 7538) — ISM: jets and outflows — Hn 
regions — stars: formation — stars: early-type — techniques: interferometry 



1. Introduction 



NGC 7538 is a large, optically visible Hn region in the Perseus arm. About 2 arcmin (1.5 
pc) southeast to the center of the Hn region is a site of active massive star formation know n 
for three radio and infrared (IR) sou rces IRS 1-3 (IMartin 19731 : IWynn- Williams et al. 1974 ). 
Located at a distance of 2.65 kpc ( iMoscadelli et al. 20091 ). the IRS 1-3 region has a total 
luminosity of 1.4 x 10 5 L Q , of which the luminosity of IRS 1 is about 8 x 10 4 L^, equiv- 
alent to a ZAMS spect ral type of 07.5 ( 1 Werner et al. 1979c ICampbell fc Thompson 1984 : 
Akabane &: Kuno 20051 ). IRS 2 and IRS 3 are two nearby sources less obscured than IRS 
1. IRS 2 has a luminosity of 5 x 10 4 L & and is associated with a compact Hn region; 
IRS 3 has a luminosity of 6 x 10 3 and is associated with an ultra compact Hn region 
f lCampbeh 1 fc Thompson 1984 : bloomer et al. 19981 : ISandell et al. 2009^ . 



IRS 1 has attracted numerous studies. Radio continuum observations show that the 
ionized gas has a do uble- lobed structure within the cen t ral 0."2 region and more extended 
halos at a 1" scale ( Campbell 1984: iGaume et al. 1995 : Sandell et al. 2009 ). and exhibits 
time variation in emission flux ( iFranco-Hernandez fc Rodriguez 20041 ) . Radio recombina- 
tion lines (RRLs) observed at cm and mm waveleng ths show extremely broad line widths , 
suggestive of expanding motion of the ionized gas ( IGaume et al. 1995c iKeto et al. 20081 ). 
In high-angular-resolution mid-IR images, a shell-like structure extending from IRS 1 to 
about 4" to the northwest is interpreted as dust emission on the wall of the outflow cavity 
( IDe Buizer fc Minier 20051 ). i^'-band images obtained with speckle interferometry show a 
fan-shaped nebula extending from IRS 1 toward the northwes t, and the emission is predom- 
inantly attributed to scattered light from an outflow cavity (IKraus et al. 20061 ). On larger 
scales, single-dish CO observations reveal a b ipolar outflow with a northwest-southea st (NW- 
SE) orientation which is centered at IRS 1 (IKameya et al. 1 989: D avis et al. 19981) . IRS 1 
is associated with various types of molecular masers (e.g.. iPratap et al. 19891 . and references 
therein). Pestalozzi et al. (2004) modeled a group of methanol masers detected toward 
IRS 1 as tracing an edge-on Keplerian disk around a central st ar of 30 Mm, although the 



orientation of a possible disk around IRS 1 is controversial (e.g., IDe Buizer &: Minier 2005 
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Klaassen et al. 20091 ). 

Existing high-angular-resolution observations toward NGC 7538 IRS 1 mostly focus 
on hot dust and ionized gas in the close vicinity of the young O star, while kinematics 
and physical conditions of molecular gas in the dusty envelope is less clear. Moreover, 
the molecular outflow in this region has not been observed at high angular resolutions; a 
pioneering interferometric study performed by Scoville et al. (1986) is of moderate resolution 
and image quality. In this paper we present a 1.3 mm waveband study of the NGC 7538 
IRS 1 region using the Submillimeter Arra}0, aimed at understanding physical conditions 
and kinematics of the gas and dust cocoon surrounding the young O star and the nature 
of the molecular outflow. In its course we discovered clustered star-forming cores deeply 
embedded within the parent cloud of IRS 1-3 and multiple molecular outflows. A description 
of observations and data reduction is presented in Section [2J Observational results are 
presented in Section [3j followed by discussions on their implications in Section HI Finally, a 
summary of the main findings is given in Section [5j 



2. Observations and Data Reduction 

2.1. SMA Observations 

The SMA observations were undertaken on 2007 July 13 with eight antennas in the 
Compact-north configuration. The weather was very good during the observation, with 1.5 
to 2.5 mm precipitable water vapor (PWV), corresponding to atmospheric optical depth at 
225 GHz, r 2 25GHz, around 0.08 to 0.12. To cover the entire outflow we observed two fields 
centered at (R.A., decL) J20 oo = (23 h 13 m 43. s 7 5, 61°28'21"5) and (23 h 13 m 46. s 75, 61°27'59."5), 



guided by previous single-dish observations (IDavis et al. 19981 ) . The 2x2 GHz correlator 



was configured to cover rest frequencies about 219.2 to 221.2 GHz in the lower sideband 
and about 229.2 to 231.2 GHz in the upper sideband, with a uniform spectral resolution 
of 812.5 kHz (~1.1 kms -1 ). We observed 3c273 and Callisto for bandpass calibration and 
MWC349, for which a flux of 1.65 Jy at 225 GHz was adopted, for flux calibration. Time 
dependent gains were monitored by interleavingly observing J0102+584 and BL Lac every 30 
mins. The primary calibrations were carried out with the IDL MIR package, which allows 
to determine antenna gains (versus frequency or time) by simultaneously fitting multiple 
calibrators. The calibrated visibilities were exported to MIRIAD for further processing. A 



lr rhe Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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continuum database was constructed from line-free channels using the MIRIAD task UVLIN. 
We self-calibrated the continuum solving for phase corrections with decreasing time intervals 
(5, 2, 1 min), and then solving for phase and amplitude corrections for the last iteration. 
The solutions from the self-calibration were applied to spectral line data as well. 

In Figure [TJ the predominant core MM1 is more than 10 times brighter than MM2-9. 
As the dynamical range of an interferometric image is limited by incomplete (u, v ) sampling, 
we assess the reliability of MM2-9. We imaged JO 102+584 and BL Lac and found that each 
of the two quasars appears as a point source at the phase center, and the position agrees with 
the nominal coordinate to < 0.1". We also found that features associated with MM2-9 are 
identifiable in the continuum and spectral line images without self-calibration. Furthermore, 
we obtained a model of MM1 by restricting CLEAN (a deconvolution task in MIRIAD) to 
a polygon closely around MM1, subtracted "MM1" from the calibrated visibilities using the 
task UVMODEL, and imaged the residual visibilities; MM1 was completely invisible while 
MM2-9 were clearly seen in the image, indicating that none of the fainter cores was falsely 
generated due to sidelobes of MM1 during deconvolution. 

The final maps have a synthesized beam of ~3" x 2", and r.m.s. sensitivities of ~5 mJy 
for continuum and ~35 mJy per 1.2 kms -1 for lines. 



2.2. Short Spacings in CO and 13 CO (2-1) 

In an interferometric observation, structures more extended than 1.2A/6 m i n , where A is 
the observing wavelength and 6 min the shortest length of the projected baselines, cannot be 
effectively sampled. This filtering effect is significant for CO and 13 CO maps, particularly at 
velocities around the cloud veloci ty. Adding in single-dish observations is a way to overcome 



this effect (e.g., lOiu et al. 20091 ). For CO (2-1), we combined the SMA data with data 
obtained from the James Clerk Maxwell Telescope (JCMT) by Davis et al. (1998). For 
13 CO (2-1), we carried out on-the-fly mapping observations with the Caltech Submillimeter 
Observatory^ (CSO) 10.4 m telescope on 2008 June 11. During the last observation, the 
weather was reasonably good with about 4 mm PWV and T225GHZ around 0.2. We obtained 
a 13 x 13 grid map centered on (R.A., decl.) J20 oo = (23 h 13 m 43. s 8, 61°28'22."0). The grid cell 
spacing is 10", roughly 1/3 of the FWHM beam at this frequency. The effective integration 
time toward each cell is about 10 s. The AOS spectrometer has 1024 channels throughout the 
50 MHz bandwidth, thus a spectral resolution of 48.8 kHz (or 0.066 kms -1 ). The data were 
reduced using the standard CLASS package. The combination of the SMA and single-dish 



2 The Caltech Submillimeter Observatory is supported by the NSF grant AST-0229008 
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data was conducted in MIRIAD following a procedure outlined in Zhang et al. (1995). The 
final CO and 13 CO maps were smoothed to 2 kms" 1 per channel. 



Results 



3.1. 1.3 mm Continuum Emission 



Figure [T] presents the 1.3 mm continuum map, which reveals a cluster of nine mm 
cores, namely MM1 to MM9 in order of peak fluxes. Measured parameters of these cores 
are listed in Table [TJ The brightest core MM1 coincides with the luminous IR and radio 
source IRS 1. The 1.3 mm continuum flux has a contribution of ~ 1.5 Jy from the free-free 



emission extrapolated from cm measurements (ISandell et al. 20091 ). and the dust emission 
flux is about 2.1 Jy. The 1.3 mm emission from MM2 to MM9 is entirely dominated by the 
dust emission, since none of them is seen in cm continuum. In contrast to IRS 1, neither 
IRS 2 nor IRS 3 appears to be embedded within a dense dusty core. 

Assuming thermal equilibrium between the dense gas and dust, we adopt a temperature 
of 245 K for MM1 derived from a local thermodynamical equilibrium (LTE) model of the 
CH 3 CN emission ( Section 13.3. ip. Cores MM2-9 have much lower temperatures; we estimate 
temperatures from NH3 data (jZheng et al. 200ll ) for MM6-9 (Section l3.3.2l) and adopt an av- 
eraged temperature of 40 K for MM2-5 (Section 14.31) . W e then estimate gas masses o f all the 
cores by adopting the dust opacity K225GHZ = 1 cm 2 g -1 (jOssenkopf fc Henning 1994T). whic h 



is equivalent to an opacity index j3 of 1.375 for k v = 10 ( 1 2 th z ) /3 cm2 § _1 Mildebrand 19831 ). 
and a canonical gas-to-dust mass ratio of 100. The derived core masses are listed in Tabled] 
The error bars of the core masses presented in Table [1] account for errors in measuring the 
continuum fluxes and calculating the temperatures. For MM1, the uncertainty in extrapo- 
lating the free- free emission is also taken into account. For MM2-5, an averaged uncertainty 
of 35% is assumed for the adopted temp erature. There is a sy stemic uncertainty, 9%, in- 
troduced by the distance error (0.12 kpc. iMoscadelli et al. 20091 ). Other uncertainties come 
from possible variations of the gas-to-dust mass ratio and the dust opacity; the latter can 
be constrained with future high angular resolution observations at different wavelengths. 



3.2. Molecular Outflows 

Figure |2] shows the velocity integrated CO (2-1) emission obtained from the combined 
SMA and JCMT data sets. Velocity channel maps of the emission are shown in Figure |3j 
Overall the CO emission shows spatially separated blue- and redshfited lobes in a NW-SE 
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orientation, whereas the detailed structures are ascribed to multiple outflows rather than a 
simple bipolar outflow seen in previous single-dish observations. 

In the integrated map, the blueshifted lobe appears as a curved filament originating 
from MM1 and diverging into three "branches" toward the northwest, north, and northeast. 
However, from the orientation, the northeastern branch with a position angle (PA) of about 
40° is more likely arising from MM4 or IRS 3; from the channel maps at —78 to —64 
kms -1 , the northern branch with a PA of about 0° may originate from IRS 3 or MM4 as 
well. In addition, at least four more dusty cores, MM2, MM6, MM7, and MM8 appear to 
be associated with blueshifted CO outflows: an elongated feature extending from MM6 to 
the southeast is seen in the integrated map and blueshifted channel maps of >— 78 kms -1 ; 
emission around MM8 and extending to the northwest is seen in the integrated map and 
blueshifted channel maps of >— 74 kms -1 ; emission associated with MM2 is seen as a knot 
to its immediate southwest at —78 to —72 kms -1 and as a compact elongation at —70 to 
—68 kms -1 ; compact and faint emission extending from MM7 to the southwest is visible at 
—70 to —66 kms -1 . 

The redshifted emission appears as a large-scale "arc" in the integrated map, while the 
detailed structures can be attributed to at least two cores, MM1 and MM7. MM1 is pre- 
sumably the powering source of the emission extending from around MM1 to a bright clump 
about 20" (0.26 pc) to the southeast; the clump is encompassed by a C-shaped structure 
seen in C 18 (2-1). MM7 is likely the central source of the CO emission in its close vicinity: 
a knot about 7" (0.09 pc) to the northeast seen at —38 and —36 kms -1 is most likely the 
redshifted counterpart of the compact feature seen at —70 to —66 kms -1 ; a southwestern 
clump identifiable in the integrated map and channel maps of >— 48 kms -1 is probably 
associated with MM7 as well. 

Figures S] and show the velocity integrated and channel maps, respectively, of the 
13 CO (2-1) emission, obtained from the combined SMA and CSO data. In general the 13 CO 
emission shows a morphology similar to that of the CO emission, but within a narrower 
velocity range. The MM1 NW-SE outflow appears more predominant in the 13 CO emission. 
Although fainter, the PA~ 40° and 0° branches of the blueshifed emission are both identifi- 
able. The redshifted emission of the MM7 outflow is clearly seen; in particular in the channel 
map of —48 kms -1 it shows a bi-conical structure. 

A schematic view of the multiple outflows detected here is shown in Figure [6j which 
presents a tentative anatomy of the complicated outflow structures seen in CO and 13 CO 
(2-1). A more detailed discussion on the MM1 outflow is presented in Section |4"72"1 



Following Qiu et al. (2009), the outflow parameters, including the mass (M out ), mo- 
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mentum (Pout), energy (E out ), dynamical timescale (tdyn), mass outflow rate (M out ) and 
momentum outflow rate (Pout), are calculated. The optical depth effect as a function of 
velocity and spatial distribution is corrected by comp aring; the CO with 13 CO fluxes and 
adopting a C-to- 13 C ratio of 79 (IWilson &: Rood 19941 ) . The derived parameters are listed 
in Table |2j To avoid contamination from ambient diffuse gas, we conservatively calculate 
the mass of the outflowing gas at velocities of <— 64 and >— 48 kms -1 , which are mor e 
than 6 kms -1 away from the systemic velocity of —57.4 kms -1 (Ivan der Tak et al. 20001 ) . 
However, from Figure El the emission at —62 and —50 kms" 1 has significant contribution 
from outflows. If we choose the lowest outflow velocities to be —62 and —50 kms -1 , the gas 
mass amounts to 110 M & . Hence, the total outflow mass should be in the range of 50 to 100 
M . 



3.3. High-density Tracing Molecular Lines 

3.3.1. Emission and absorption toward MM1 



MM1 shows mo lecular chemistry typical of a hot molecular core (HMC, ICesaroni 2005 



Beuther et al. 20091 ) . Excluding CO and 13 CO (2-1), around 70 lines arising from 15 molec- 
ular species are detected toward MM1 with peak fluxes above 6a (Table [3]). Spatially unre- 
solved emission from complex O-bearing molecules, e.g., methyl formate (CH3OCHO) and 
ethanol (C 2 H 5 OH), and N-bearing molecules, e.g., methyl cyanide (CH 3 CN) and isocyanic 
acid (HNCO), is clearly seen. In particular, the number of detected lines is largely dominated 
by O-bearing molecules, e.g., CH 3 OH, CH 3 OCHO, and C 2 H 5 OH. 



The CH3C N emission is an excellent temperature and density probe (lAraya et al. 2005 



Cesaroni 20051 ). The CH 3 CN lines detected in MM1 cover a wide range of upper level en- 
ergies above the ground state, P up (69-526 K). Following Qiu & Zhang (2009), we perform 
a simultaneous fitting to all the detected K components of the CH 3 CN (12-11) emission, 
where the temperature, column density, and source size are determined by means of grid 
search %-square minimization. In Figure [TJ the best fit model agrees very well with the 
observation, and yields a temperature of 245^25 K, a column density of 5±Jxl0 16 cm -2 , 
and a source size of 1 450^0 AU. As the fitting is highly nonlinear, uncertainties are deter- 
mined from a Monte Carlo analysis accounting for the la noise and the ~10% uncertainty 
of flux calibration. Compared to the rotation diagram method, the fitting impro ves the 



parameter derivation by solving for the optical depths ( iGoldsmith fc Langer 19991 ). which 
range from 0.16 for the K=8 line to 2.99 for the K=3 line in the best fit model. The fit- 
ting, however, simply assumes uniform gas. Density and temperat ure gradients certainly 



exist in a gas and dust envelope surrounding a (proto)star (e.g., IScoville fc Kwan 1976 
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Osorio et al. 20091 ). The best fit model thus represents a characteristic estimate. The beam- 
averaged CH 3 CN column density, 3xl0 15 cm -2 , can be compared to the H 2 column density 
(1.4 xlO 24 cm -2 , derived from the dust continuum peak) for an estimate of the CH 3 CN 
abundance, which is 2.1xl0~ 9 , in agreement with typical values of 10~ 10 -10~ 8 measured in 



other HMCs (e.g., lHatchell et al. 19981 : IChen et al. 2009 : IZhang et al. 20071 ). Finally, the 



C2H5OH (131,13-120,12) (i/ =220601.93 MHz) line in Figure 7 is often identified as CH 3 13 CN 
(12 3 -11 3 ) (z/ =220599.98 MHz) in the literature, which would imply a CH 3 CN (12 3 -11 3 ) 
optical depth of ~42, abnormally high for such a high-density and warm gas tracer. If the 
optical depths were that high the low K components of the CH3CN (12-11) emission would 
be entirely optically thick and show flattened line peaks; this clearly disagrees with the 
observed spectra. Here we believe the emission arising from C2H5OH (13i,i 3 -12o,i2). 

A secondary emission feature around —50 kms" 1 is seen in Nitrogen- and Sulfur-bearing 
species, i.e., CH 3 CN, HNCO, OCS, SO, and S0 2 . The feature is not detected (< 6cr) in lines 
of H^CO, CH3OH, or any other large O-bearing molecules. The emission is probably arising 
from an interaction between the outflowing material and a clump residing in the far side of 
the central young O star. The lines in Nitrogen- and Sulfur-bearing species could be related 
to shock-induced chemistry, while the deficiency in lines of CH 3 OH, H^CO, and larger O- 
bearing molecules could be d ue to an initially low abundance of the "key" species CH 3 OH 



and H 2 CO for the clump (e.g.. ICharnley et al. 19921 ). The existing observations cannot shed 
further light on this issue. 

Thought to be a probe of inward motion of dense gas surrounding a forming star, molec- 
ular lines with redshifted absorption observed at high-angular-resolutions have been rep orted 



in a very few high-mass star- forming cores (e.g.. lZhang et al. 19971 ; iBeltran et al. 20061 ). To- 
ward MM1, an inverse P-Cygni profile with redshifted absorption (<— 6cr) is seen in lines 
of C 18 (2-1), SO (5 6 -4 5 ), HNCO (10 , 10 -9o, 9 ), CH 3 OH (8_ lj8 -7 , 7)E and (8 , 8 -7 1; 6)E, 
and CH 3 CN (12-11) K=2,3 (see, e.g., Figure [HD- In these lines, the absorption comes from 
dense gas lying on the near side of the central continuum source, absorbing bright contin- 
uum emission and moving toward the source (away from the observer). Meanwhile, the 
emission is attributed to dense gas residing on the far side of the central source and moving 
toward the observer; within a 3" x 2" beam, dense gas with velocities largely along the plane 
of sky (i.e., orthogonal to the line of sight) can contribute to the emission as well. The 
lines showing redshifted absorption have E up spanning a range of 16-133 K. Roughly speak- 
ing, the absorption in lines of higher excitation (e.g., the CH3OH and CH3CN lines) may 
trace inner, warmer gas with velocities mostly along the line of sight, hence appears more 
biased toward redshifted velocities. In addition, analytical models of a dense core under 
self-similar gravitational collapse shows that the infall velocity scales with the radius as r 



-0.5 



(ILarson 19721 ; IShu 19771 ). which may partly contribute to the redward shift of the absorption 
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in lines of higher E up . Assuming the collapse proceeds in a spherical geometry, the accretion 
rate can be estimated following M acc = 47ir 2 pV in f a u, where M acc is the mass accretion rate, 
r the radius of the infalling envelope, p the mass density, and Vi n f a u the infall velocity. From 
the dust continuum peak the H 2 number density is about 10 7 cm -3 averaged over r ~ 3000 
AU. The measured infall velocity is about 2.5 kins" 1 . As an order-of-magnitude estimate, 
M acc amounts to 3 x 10~ 3 M Q yr _1 . However, the accretion may considerably deviate from 
a spherical geometry, in particular for the inner part closely surrounding the central source. 
The free-free emission in IRS 1 shows a double-lobed morphology at the center, suggest- 
ing that the accretion proceeds along a flattened structure perpendicular to the ionized gas 
elongation. If the collapse occurs inside a solid angle il, M acc scales to ,0/(471") • 3 x 10~ 3 
Moyr- 1 . 

3.3.2. Emission from new cores MM2 to MM9 

In addition to MM1, the dust continuum map reveals eight new cores, MM2-9. The 
C 18 (2-1), SO (56-4 5 ), and CH3OH (8-i^-Toj) E lines show emission from part of these 
new cores and extended structures, and are not significantly affected by outflows or ambient 
diffuse gas. Figure [9] shows moment-0 (velocity averaged) maps of these lines. 

The C 18 (2-1) emission is dominated by two dense clumps residing on both sides of 
an absorption feature around MM1: the western clump peaks at MM2 and shows extension 
from MM5 to MM4; the eastern clump encompasses MM3 and MM6. While MM2 coincides 
with the brightest C 18 peak, cores MM3-6 are more or less offset from local C 18 peaks. 
Another interesting feature in C 18 is a C-shaped shell seen in the southeast of the mm 
cores. In comparison with the CO and 13 CO (2-1) emission, the C 18 shell lies immediately 
ahead of the brightest and furthest CO/ 13 CO clump in the redshifted outflow (see Figure [2] 
and Figure H]). This suggests that the C-shaped structure forms from ambient gas which is 
compressed by the outflow from MM1. 

Apart from emission and absorption toward MM1, the SO (5g-4s) map shows bright 
emission around MM2, MM3, MM6, MM8, MM9, as well as faint extensions toward MM4 
and MM5. The SO map also shows clumps without a counterpart in the dust continuum; for 
example, two bright clumps to the southeast of MM6 are not seen in dust continuum or any 
other line. It is unclear whether the SO emission in the two clumps is associated with the 
MM6 outflow or arising from protostellar cores whose dust continuum is below the detection 
limit of our SMA observation. 

Compared to C 18 (2-1) and SO (5 6 -4 5 ), the CH 3 OH (8_i )8 -7 , 7 ) Eline has a consider- 
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ably higher E up (89 K) and appears to trace more compact structures. In addition to MM1, 
cores MM3, MM6, MM7, MM8, and MM9 are clearly detected in this line. The CH 3 OH 
peaks approximately coincide with the dust continuum peaks, indicating the presence of 
internal heating within each of the detected cores. Relatively faint emission associated with 
MM5 is detected as well. 

A serendipitous discovery from the CH 3 OH (8-1^—70,7) E map is a new 229.8 GHz 
maser at (R.A., decl.)j 2 ooo = (23 h 13 m 45. s 52, 61 27'41."5), which is seen as unresolved and 
bright emission with a peak flux of ~1.64 Jybeam -1 and a narrow line width of less than 2 
kms" 1 (Figure |9b and Figure [TO]) . It is difficult to assess whether this maser spot is associated 
with any of the mm continuum sources which are more than 0.2 pc away. Another possibility 
is that the maser is associated with a compact CH3OH (8_i i8 -7 0i 7) E and SO (5 6 -4 5 ) clump 
about 5" to the north. The clump probably traces a low- to intermediate-mass protostellar 
core whose 1.3 mm dust emission is below the detection limit of our SMA observation. 
More interestingly, the maser spot lies right at the tip of a compact CO outflow (Figure |2]), 
suggesting that the masing line is pumped by the outflow shocks, which is consistent with 
the class I property of the t ransition. There has been very few r eport on the detection of 



229.8 GHz methanol maser (jSlysh et al. 2002t lOiu fc Zhang 20091 ). probably because of its 



high frequency hence not covered by previous methanol maser surveys. 

In the VLA NH 3 (1,1) and (2,2) data taken by Zheng et al. (2001), emission in MM6- 
9 is clearly seen in both transitions (see Figure [TT]) . while an investigation of emission in 
MM2-5 is severely affected by predominant absorption arising from MM1. As there is no 
clear detection of a satellite line in Figure [TTJ we assume the emission to be optically thin, 
and deduce rotational temperatures, Tr, from the intensity ratio of the two transitions fol- 
lowing T R = — 41.7/in(0.282T 2j 2; m /T 1>1;m ), where T 11;m and T 2 ^-m are the brightness temper- 
atures of the main h yperfine components of the NH 3 (1,1) and (2,2) transitions, respectively 



( iHo fc Townes 19831 ); the formula assumes equal line width for the two transitions, which 
is probably a good approximation but difficult to verify given the moderate sensitivity and 
spectral resolution. We obtain T R of 29 ± 13, 34 ± 12, 33 ± 14, and 27 ± 7 K for MM6-9, 
respectively; the un certainties account for the la noise of the spectra. Correcting for the 



depopulation effect (IDanby et al. 19881 ). we estimate the kinetic temperatures from T R to 



be 43, 58, 54, and 38 K for MM6-9, respectively. 
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Discussion 



4.1. Nature of the MM1 core 

The MM1 core has a size of about 0.05 pc, an H 2 number density of order 10 7 cm -3 , 
a characteristic temperature of 245 K, and rich chemistry in complex organic molecules. 
All these indicate that MM1 is a typical HMC On the other hand, the presence of bright 
IR and radio emission suggests t hat the core is at a relatively more evolved stage com 



pared to many oth er HMCs (e.g., IChen et al. 2009 : IZhang et al. 20071 iGirart et al. 2009 



Qiu fc Zhang 20091 ). The core also shows particularly rich chemistry in large O-bearing 
molecules (e.g., CH3OCHO and C2H5OH), consistent with a relatively more evolved stage. 

According to stellar structure models of Schaller et al. (1992), the central young star 
embedded within MM1 has a mass of > 25 M . Redshifted absorption seen in several 
lines toward MM1 implies the presence of an accretion flow onto the central star or star- 
disk system which may further increase the mass of the young O star. The ~20 M & core, 
however, is not equivalent to the gas reservoir of the accretion. On the one hand, My- 
ers (2009) found that clustered and high-density regions form high-mass protostars from 
both core and environment gas. There are also claims that t he accretion onto a formin g 
O star may be fed by global infall of the parent cloud (e.g., iGalvan-Madrid et al. 20091 ). 
For NGC 7538 IRS 1, it is unclear whether there is a global infall pertaining to the accre- 
tion. Previous NH3 (1,1), (2,2) observations revea led blueshifted absorption toward IRS 1 



( IWilson et al. 19831 ; iHenkel et al. 1984 ; iKeto 199ll also see Figure [H]) , which is attributed 
to expanding or outflowing gas moving toward the observer. The NH 3 absorption is opti- 
cally thick (r ^ 1, Henkel et al. 1984; Keto 1991), presumably tracing gas mostly residing 
in outer layers surrounding the HMC. Thus the bulk of the outer gas appears to be in an 
outward motion. On the other hand, various mechanisms, e.g., outflows, heating and ioniza- 
tion, and accr etion by nea rby protostars, can disperse dense gas and conspire to terminate 
the accretion ( iMyers 20091 ). In particular, Keto (2007) described a three-stage evolutionary 
sequence of an Hn region around a massive young star. In his model, the dynamics and 
morphology of the ionized gas are determined by the ratio of the ionization radius, Ri, and 
the gravitational radius, Rb, defined by GM Jf /2c 2 , where M* is the mass of the star and c 
is the sound speed. With M* = 25 M and c = 10 kms -1 (for ionized gas temperature 
of 1 4 K), Rh is 110 AU. The innermos t ionized gas around IRS 1 has a radius of ~500 
AU JCampbell 19"83 baume et al. 19951 : ISandell et al. 2009f ). Thus, for IRS 1, Ri > R b , 
corresponding to the second to third stages in Keto's model, that is, the ionized gas moves 
outward and progressively expands toward the equatorial plane. In this picture the accretion 
is confined to a narrow range of angle around the disk and is close to termination. 
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4.2. The MM1/IRS 1 outflow 



Multiple outflows are seen in the CO and 13 CO maps. However, the bulk outflowing gas 
is centered at the brightest core MM1 which harbors the IRS 1 young O star. As outlined 
in Figure O the blueshifted lobe of the MM1 outflow is dominated by a slightly curved 
filament with a PA of about —50°. In sub-arcsec resolution observations, an elongated feature 
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be traced back to about 0."2 (~500 AU) from the central source; within r 
ionized gas shows a double-lobed morphology in a north-south orientation. In addition, 
near-IR observations toward IRS 1 reveals a fa n-shaped structure (opening angle ~ 90°) 
extending from IRS 1 to the northwest (~ 10", iKraus et al. 20061 ). On a larger scale, the 



CO/ 13 CO filament is seen projected against a dense clump revealed by the NH3 emission 
(see Figure E]). The most likely scenario capable of incorporating the inner ionized gas with 
a north-south orientation and the northwestern structure extending from 500 AU to 0.5 pc 
from the central source seems to be a wide-angle wind from IRS 1 which carves an outflow 
cavity. The axis of the wide-angle wind is mostly along the north-south orientation and 
slightly inclined to the west. The northwestern structure observed in radio, mid-IR, and 
in CO and 13 CO (2-1) emission traces the western wall of the outflow cavity. The lack of 
appreciable emission from the e astern wall is likel y ascribed to the deficiency of dense gas 
to the northeast of IRS 1 (e.g., IZheng et al. 200ll ). Indeed, from a close inspection of the 



ionized gas and hot dust immediately around IRS 1, faint extension to the northeast does 
exist as well (Figure 2 in Campbell 1984; Figure lb in Gaume et al. 1995; Figure 2b in 
De Buizer & Minier 2005). The CO/ 13 CO structure originating from the vicinity of IRS 3 
and MM4 and extending to the north (PA~0°) is puzzling; in Figure |2] and Figure HJ the 
structure seems to be part of the MM1 outflow but with an orientation very different from 
the inner part (i.e., the P A of the outflow yarying from —50° to 0°). In dust continuum 



or NH 3 maps of the re gion (IZheng et al. 200 ll ; ISandell & Sievers 20041 : iReid fc Wilson 2005 



Pestalozzi et al. 20061 ). there appears to be no dense gas immediately to the west of the 
structure, excluding the possibility of being confined by ambient dense gas. From Figure |6l 
the wide-angle wind seems to be expanding into an ambient dense clump seen in NH3 and 
as a consequence of this interaction, part of the gas is deflected to the north, forming the 
PA~0° structure. However, with the existing data it cannot be ruled out that the structure 
is associated with IRS 3 or MM4. 

It is not straightforward either to unambiguously identify and interpret the redshifted 
lobe of the MM1 outflow. In the wide-angle wind scenario proposed above, the structure 
mostly consists of ambient gas being entrained or swept up by the wide-angle wind and 
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tracing the eastern wall of the southern outflow cavity. The slightly distorted appearance 
of the structure could be partly caused by interaction with other outflows (e.g., the MM7 
outflow). An interesting feature is the C-shaped shell seen in C 18 0. It lies ahead of the 
furthest and brightest CO/ 13 CO clump and features dense gas being compressed by part of 
the leading front of the wide-angle wind. It is puzzling that the western wall of the southern 
outflow cavity is unseen in CO or 13 CO, as at least in projection dense gas to the southwest 
of MM1 is seen in NH 3 and dust continuum. One possibility is that in a three-dimensional 
picture the southern lobe of the wide-angle wind is expanding into a medium without dense 
gas in the west and cannot produce appreciable emission. 

Hollenbach et al. (1994) studied the photoevaporation of disks around massive young 
stars and developed a model of the photoevaporative wind from the disk. Lugo et al. (2004) 
performed a parametric investigation of such a wind in NGC 7538 IRS 1 by modeling the free- 
free emission as arising from the wind and found a mass- loss rate of 1.1 x 10~ 5 Mp yr -1 . As the 



speed of the photoevaporative wind is ~10-50 kms -1 (jHollenbach et al. 1994 ). comparable 
to that of CO outflows, the mass loss rate of the wind seems too low to drive the NGC 7538 
IRS 1 outflow, whose mass loss rate is of order 10 -3 M & yr -1 (Table[2]). The photoevaporative 
wind may play a role in driving the IRS 1 molecular outflow, however, the large amount of 
gas mass, momentum, and energy in outflowing gas requires a more energetic wind as its 
driving engine. In a momentum-driven picture (working like a snow-plow), a fast, ~1000 
kms -1 wind with a mass loss rate of 10" 5 M yr 1 is able to drive the observed molecular 
outflow; in an energy- driven situation (working by ram pressure), the required mass loss 
rate is much lower, i.e., 10~ 7 M yr -1 , provided the wind speed reaches 1000 kms -1 . Stellar 
winds from O stars and possible accretion-driven winds (e.g., a high- mass analogue of X- 
winds in low-mass protostars) are expected to have velocities reaching ~1000 kms -1 , and 
energy-wise, are likely to drive the massive and energetic molecular outflow. 



4.3. A massive cluster in the making 



In single-dish submm and mm continuum maps of the NGC 7538 com plex, the most re- 
markable cloud has a size of about 1 pc and s lightly extends to the southeast (ISandell fc Sievers 20041 ; 
Reid fc Wilson 20051 ; iPestalozzi et al. 20061 ) . Star formation activity within this cloud has 
long been thought to be concentrated around IRS 1-3, and IRS 1 is the source that attracted 
most of the interest. Our SMA 1.3 mm continuum observation uncovers a total of nine cores 
with a spatial distribution roughly following the shape of the cloud. Apart from MM1, which 
harbors the IRS 1 young O star, cores MM2-9 are new detections. The detection of these 
dusty cores is interesting in the context of cloud fragmentation and cluster formation. 
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For MM1, with a number density of 1.7 x 10 7 cm -3 (averaged within a beam and 
assuming a spherical geometry) and a temperature of 245 K, one derives a thermal Jeans 
mass (Mj) of 17 M , which is comparable to the measured core mass and compatible with 
the scenario that heati ng from the forming star may have helped t o suppress fragmentation 



of the envelope (e.g., iKrumholz 2006t iKrumholz fc McKee 20081 ) . However, all the new 



cores are very young, apparently lacking significant heating. These cores are not detected 
in CH 3 CN, and only part of them are seen in one transition of CH 3 OH (see Figure [9]). For 
MM6-9, the NH 3 data reveal temperatures of order 50 K (Section GE£2J. For MM2-5, a 
rough upper limit may be obtained based on the absence of CH3CN emission, which is a 
typical thermometer of HMCs. For example, cores MM2-5 have masses close to MM1; if the 
abundance ratio and spatial distribution of CH3CN in MM2-5 resembles that in MM1, at 
a temperature of 100 K, the brightness temperature (averaged within a beam of 3" x 2") of 
the CH 3 CN (12-11) if =0-3 emission in MM2-5 would be about 2-4 K, which could have 
been detected at a > 15<r level. Alternatively, with the CH3CN (12-11) if-ladder observed 
with the SMA, Qiu & Zhang (2009) derived a temperature of 110 K for a HMC in the HH 
80-81 massive star-forming region; if the HH 80-81 HMC were at a distance of 2.65 kpc 
and observed with a 3" x 2" beam, the brightness temperature of the K=0,1 components 
would have been around 1 K; cores MM2-5 in NGC 7538 have masses comparable to that 
of the HH 80-81 HMC, while the CH 3 CN emission is not detected (a ~ 0.13 K). All this 
indicates that MM2-5 have temperatures well below 100 K. Sandell & Sievers (2004) found 
a dust temperature of 40 K for sub mm condensations around IRS 1. This probably provides 
a rough estimate of the averaged temperature of MM2-9, as part of them roughly coincide 
with the submm condensations. 

MM2-9 have number densities of ~ 0.4-1.2 x 10 7 cm -3 , or mass column densities > 1 
gem -2 , matching a proposed threshold of 1 g cm" 2 for molecular clou ds capable of avoiding 



fragmentation and forming massive stars ( IKrumholz fc McKee 20081 ). But the masses of 



the new cores are well above Mj, implying the importance of supersonic turbulence and/or 
magnetic field in c loud fragmentation and possible star formation within the new cores (e.g., 



Zhang et al. 20091 ). While an investigation of the magnetic field awaits new observations, a 
glimpse of turbulence can be obtained by looking into the emission lines arising from these 
cores. Figure DJ shows spectra in C 18 (2-1), SO (5 6 -4 5 ), and CH 3 OH (8_ 1)8 -7 0j7 ) Eat 
the peaks of MM2-9. A gaussian fitting is performed on lines with sufficient signal-to- noise 
ratios. The derived line widths of ~2.5-4 kins" 1 are more than 5 times larger than the 
thermal line width. This strongly suggests the presence of supersonic turbulence which may 
play an important role in suppressing fragmentation, if the line widths are not dominated 
by kinematics related to star formation processes (e.g., outflow, infall, and rotation). 

Information on potential ongoing star formation within the new cores can be obtained 
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from their associations with CO outflows, H2O masers, and warm and dense gas conden- 
sations. Cores MM2, MM3, MM4, and MM7 are all associated with H 2 maser emission 
(Figured]), indicative of active star formation. Cores MM2, MM6, MM7, MM8, and prob- 
ably MM4 as well, are associated with CO outflows (Figure [6] and Section l3T2~j) . implying 
the presence of protostars within these cores. MM9 shows the faintest dust continuum; al- 
though lacking apparent association with an outflow or H2O maser spot, it is very likely 
internally heated by an intermediate-mass protostar as it shows centrally peaked emis- 
sion in CH3OH and SO (Figure E]). There is no clear signature of star formation toward 
MM5, which also shows a relatively flattened profile in dust emission. It is more likely 
a gas and dust clump exte rnally heated by IRS 2 (probably through a stellar-wind bow 
shock, iBloomer et al. 19981 ). rather than an internally heated star-forming core. In short, 
seven out of the eight new cores appear to be forming stars. The star-forming cores have 
masses of ~ 5-24 M & and high densities, and most of them show bright, centrally peaked 
CH3OH emission. Although CH3OH lines are seen i n hot corinos (warm and dense envelopes 
around low-mass protostars, ICeccarelli et al. 20071 ). the emission could not be detected if 
the hot corinos were at d i stances > 2 kpc as compared to the c ase of IRAS 16293-2422 
( Ivan Dishoeck et al. 1995t iKuan et al. 2004t IChandler et al. 20051 ). Thus, the new cores 
seen in CH3OH are expected to be forming intermediate- to high-mass stars. MM4 is not 
detected in CH3OH and its immediate northwest is deficient of emission in C 18 and SO; 
molecular line emission from this core is probably hampered by UV radiation from IRS 
3, which has begun to ionize the surrounding and is not heavily obscured by a dusty co- 
coon. MM2 is the most massive of the new cores and coincides with the brightest C 18 
peak, but is not seen in CH 3 OH. With a relatively low E up (16 K), the C 18 (2-1) emission 
preferentially traces c olumn density enhancements rather than temperature increase (e.g., 



Wyrowski et al. 19991 ). The CH 3 OH (8_i ( g— 7o, 7) E emission (with a much higher E uv of 89 



K) is more sensitive to the heating. Therefore, MM2 is presumably the youngest of the new 
star-forming cores. 

The NGC 7538 molecular complex has been t hought to exhibit sequent i al massive star 
formation propagating from northwest to southeast (jElmegreen fc Lada 19771 ; IWerner et al. 19791 ) 
In this picture, the IRS 1-3 region lies ahead of the ionization front of the optical Hn re- 
gion and represents an intermediate evolutionary phase between the ionizing stars of the 
Hn region in the northw est and deeply embe dded massive protostars in the southeast (i.e., 
NGC 7538 S and IRS 9, IWerner et al 1979h . By zooming in the parent cloud of IRS 1-3, 
our SMA observations provide new insights into this scenario. With the newly discovered 
star-forming cores, the IRS 1 young O star is not only accompanied by more evolved sources 
IRS 2 and IRS 3, but also surrounded by a group of intermediate- to high-mass protostars 
distributed over the parent cloud. There is no pronounced propagation of star formation 
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across the cloud. In comparison with massive protostars in IRS 9 or NGC 7538 S, IRS 1-3 
are certainly more evolved, but all the newly discovered star-forming cores are at very young 
evolutionary stages. In particular MM2, which is the closest to MM1, is apparently younger 
than IRS 9, which is visible in the near- to mid-IR dWerner et al. 1979h. and eve n younger 
than NGC 7538 S, which shows bright CH 3 CN emission JSandell fc Wright 201ph . The pc- 
sized cloud appears to be in an active phase of cluster formation, with around 10 OB-type 
stars and probably orde rs of magnitude more low-mas s stars according to a typical IMF (ini- 
tial mass function, e.g., ISalpeter 1955t Kroupa 20021 ). Indeed, within the volume occupied 
by high- and intermediate-mass protostars described above, Kraus et al. (2006) in their 2.2 
//m speckle images found a number of weak, compact sources. It is interesting to speculate 
that these are already more evolved lower mass members of the same forming star cluster 
of which IRS 1-3 and the new protostars (MM2-4 and MM6-9) define the higher mass end 
of the mass distribution. Here, a comparison with the well-studied Orion Nebular Cluster 
(ONC) is instructive. As argued by Palla & Stahler (1999) and Huff & Stahler (2006), the 
highest mass ONC member stars comprising the "Trapezium" at the center of the cluster are 
~ 10 5 yr old, while the bulk of the lower mass cluster members formed 1-3 million years ago. 
In fact, one may consider IRS 1-3, MM2-4 and MM6 to be a Trapezium system in forma- 
tion. The maximal physical separation between these objects, 0.30 pc, is much larger than 
the maximal separation between members of the ONC Trapezium, 0. 047 pc, i.e., betwee n 
9 1 Ori D and 9 l Ori E, assuming a distance of 417 pc to the ONC fjMenten et al. 20071 ). 
It is, however, well within the range of other Trapezium-type systems. A Trapezium-type 
system, as defined by Ambartsumian (1954; 1955; 1958), of which 9 l Ori is the prototype, 
is a multiple star whose separations are all roughly of the same order of magnitude. While 
their cosmogonic role is somewhat unclear, it is a well-established fact that such systems are 
frequently found as the highest mass systems at the centers of clusters of lower mass stars 
(ISharpless 19541 ) . A catalog of 87 Trapezia with maximum separations, r max , within 0.2 pc 
was published by Salukvadze (1978), while Salukvadze & Dzhavakhishvili (1988) reported 15 
wider systems with r max between 1.1 and 5.7 pc. Thus the dimensions of the nascent NGC 
7538 IRS 1 Trapezium are well comparable to values generally found for such systems. 



5. Summary 

We present a 2-3" resolution study in the 1.3 mm waveband of the well-known massive 
star-forming region NGC 7538 IRS 1-3. The gas and dust core MM1, which harbors the 
IRS 1 young O star, is massive (~ 20 M ), hot (~ 245 K), and shows bright line emission in 
complex organic molecules. Redshifted absorption seen in several high-density tracing lines 
indicates that the central young O star or star-disk system is still accreting gas from MM1. 



In comparison with a theoretical model on early evolution of Hn regions (IKeto 20071 ). the 
accretion seems to be confined to a narrow range of angle about the equatorial plane. 



Neither of the two nearby sources IRS 2 and IRS 3 is associated with a dense dusty 
envelope, confirming that the two luminous sources are at most moderately obscured. 

Sensitive and high-angular-resolution observations also reveal eight new cores MM2-9 
embedded within the parent cloud of IRS 1-3. These cores have masses much larger than 
the thermal Jeans mass, indicating the importance of turbulence and/or magnetic field in 
cloud fragmentation. The cores are seen in one or more of the C 18 (2-1), SO (56-45), 
and CH3OH (8-1,8-70,7) E lines; the line widths are significantly larger than the thermal 
line width, attributed to supersonic turbulence and/or kinematics related to ongoing star 
formation. Seven out of the eight cores appear to be forming intermediate- to high-mass 
stars, of which the most massive ones appear to be members of a nascent Trapezium system. 
Characterized by the deeply embedded young O star IRS 1, two moderately obscured sources 
IRS 2 and IRS 3, and seven new intermediate- to high-mass protostellar cores, and taking 
into account presumably more unseen low-mass (proto) stars, a massive cluster is forming 
from a pc-sized cloud exposed to the ionization front of an optical Hn region. 

The CO and 13 CO maps obtained from the combined SMA and single-dish data reveals 
an outflow scenario far more complicated than a bipolar outflow seen in single-dish obser- 
vations. Multiple outflows are seen, and the MM 1 /IRS 1 outflow is dominating the mass 
of 50 M Q , momentum of 4.6 x 10 2 M Q kms -1 , and energy of 4.9 x 10 46 erg calculated for 
the gas moving at velocities of < —64 kms -1 and > —48 kms -1 . Comparing our CO, 13 CO 
maps with the existing radio and IR observations we suggest that the MM1/IRS 1 molecular 
outflow is likely driven by a wide-angle wind from the forming O star or star-disk system, 
and the wind is expanding into a medium with inhomogeneous density distribution, creating 
asymmetric and curved structures seen in CO/ 13 CO. 

We are grateful to Dr. Jingwen Wu and Dr. Ruisheng Peng for help in obtaining the 
13 CO (2-1) observations at the CSO. We acknowledge the anonymous referee for construtive 
comments which improved the paper. 
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Fig. 1. — 1.3 mm continuum emission shown in solid contours with contour levels increasing 
from 30 to 165 mJybeam -1 in steps of 15 mJybeam -1 and continuing to 3255 mJybeam -1 
in steps of 30 x (1, 2, 3, 12) mJybeam -1 ; dotted contours show negative emission with 
absolute levels the same as that of the positive. The plus symbols depict water maser spots 
from Kameya et al. (1990); three filled dots mark the positions of IRS 1-3, where IRS 1 
roughly coincides with MM1, IRS 2 is about 10" to the north and IRS 3 about 14" to the 
west. Hereafter a filled ellipse in the lower left delineates the synthesized beam at FWHM. 
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23 h 13 m 50 s 40 s 

Right Ascension (J2000) 

Fig. 2. — Integrated CO (2-1) emission obtained from the combined SMA and JCMT data. 
Blue contours show emission integrated from —78 to —64 kms -1 and the contour levels 
start from 10.2 Jy-kms -1 and increase in steps of 6.8 Jy-kms -1 ; red contours show emission 
integrated from —48 to —30 kms -1 and the contour levels start from 20.2 Jy-kms -1 and 
increase in steps of 10.1 Jy-kms -1 . The numbered squares mark dusty cores MM 1-9 (see 
Figure [1]); three thick circles denote IRS 1-3; a filled triangle marks a newly discovered 
229.8 GHz CH 3 OH maser (see Section 13.3.21 and Figures 19b. fTDl. A gray curve delineates a 
C-shaped structure seen in C 18 (2-1) (see Section [3T3T2] and Figure Ek), and dashed arrows 
outline the three "branches" of the blueshifted outflow (see Section |3~2]) . 
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61 °27'30" k ........... : . 

23 h 13 m 50 s 40 s 



Right Ascension (J2000) 

Fig. 3. — Velocity channel maps of CO (2-1) emission obtained from the combined SMA 
and JCMT data. The starting and spacing contour levels are 1.2 Jybeam -1 for channels 
from —66 to —46 kms -1 and 0.6 Jybeam -1 in outer line wings. The central velocity of each 
channel is shown in the upper left of each panel; cores MM1-9 and IR sources IRS 1-3 are 
denoted as filled squares and dots, respectively. 



-25 - 



29'00" - 



o 
o 
o 

C\J 



o 

H — ' 

75 

CD 

Q 



30" - 



6"T28'00" - 




23^13 m 50 s 



Right Ascension (J2000) 



Fig. 4. — Integrated 13 C0 (2-1) emission obtained from the combined SMA and CSO data. 
Blue contours show emission integrated from —72 to —64 kms -1 and the contour levels 
start from 2.6 Jy-kms -1 and increase in steps of 1.3 Jy-kms" 1 ; red contours show emission 
integrated from —48 to —42 kms -1 and the contour levels start from 3.6 Jy-kms -1 and 
increase in steps of 1.2 Jy-kms -1 . Symbols are the same as those in Figure [2j 
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Right Ascension (J2000) 

Fig. 5. — Velocity channel maps of 13 CO (2-1) emission obtained from the combined SMA 
and CSO data. The starting and spacing contour levels are 0.8 Jybeam -1 for channels from 
—60 to —56 kms -1 and 0.45 Jybeam -1 in outer line wings. Symbols are the same as those 
in Figure [3j 
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Fig. 6. — A schematic view of the multiple outflows seen in the CO and 13 CO emission. 
Dashed curves roughly outline the MM 1 /IRS 1 outflow; arrows approximately denote the 
orientation and extension of suggested outflows from the other cores; two dash-doted lines 
bisecting at MM7 highlight a bi-conical structure remarkable in 13 CO; a PA~0° structure is 
show in a dotted arrow for its ambiguous interpretation (see Section H~2l for details). Blue 
and red colors represent blue- and redshifted velocities of the structure, respectively. Other 
symbols are the same as those in Figure [2 The NH3 (2,2) emission is shown in filled gray 
contours (data from Zheng et al. 2001; see Section fl~2l for discussions on a dense clump to 
the west of MM4). 
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Fig. 7. — The observed spectra of CH 3 CN (12-11) toward MM1 in black histogram, overlaid 
with the best fit LTE model in thick gray. 
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-65 -60 -55 -50 -45 
Velocity (km s ) 

Fig. 8.— The spectra of C 18 (2-1), SO (5 6 -4 5 ), HNCO (10 0) io-9o,9), CH 3 OH (8_ li8 -7 0)7 ) 
E, and CH 3 CN (12-11) K=3 toward the peak of MM1. Three dashed lines denote channel 
velocities at —57, —55.8, and —54.6 kms" 1 , the peaking velocities of the absorption. The 
dotted gray line denotes the channel of —59.4 kms -1 , the peaking velocity of the emission 
in all the lines. The upper energy level of each transition is indicated in the right-hand of 
each panel. For comparison, the absorption in NH 3 (2,2) is shown in the bottom panel (data 
from Zheng et al. 2001). 
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Fig. 9. — Moment-0 (velocity averaged) maps shown in contours for the lines of (a) C 18 
(2-1), (b) SO (5 6 -4 5 ), (c) CH 3 OH (8_i i8 -7 0l7 ) E. Contour levels start from 50 mJybeam" 1 
and increase with steps of 25 mJybeam -1 for the C 18 map, start from 48 mJybeam -1 and 
increase with steps of 24 mJybeam -1 for the SO map, and start from 24 mJybeam -1 and 
increase with steps of 29 mJybeam -1 for the CH 3 OH map. Cores MM1-9 are labeled as 
filled squares and for clarity, with their nomenclature shown in the CH 3 OH map; IRS 1-3 
are denoted as three dots. 
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Fig. 10. — The velocity profile of the 229.8 GHz methanol maser discovered in Figure [9b (see 
Section [3.3.21 for details). 



- 32 - 



0.04 


MM6 


NH 3 (1,1) : 


IMM6 


NH 3 (2,2) : 


0.02 


1 _ i«l 


U ii 


_flfl nil n 


: 




MM7 

0.04 P 
0.02 ; 

S 00 ^ 

_Q -0.02 [ 

3 °- 04 : 
X 0.02 - 



■i i i- 



MM7 




i i i i i i i- 



MM8 



-0.02 ; 
0.04 - 



■i i i- 



MM8 




MM9 



0.02 [ 

0.0 
-0.02 P 



MM9 




-90 -80 -70 -60 -50 -40 -30-90 -80 -70 -60 -50 -40 -30 

Velocity (km s" 1 ) 

Fig. 11.— The NH 3 (1,1) and (2,2) spectra in MM6-9; data taken from Zheng et al. (2001) 
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Fig. 12.— The C 18 (2-1) Left), SO (5 6 -4 5 ) (Middle), and CH 3 OH (8_ lj8 -7 0i7 ) E (Right) 
spectra toward the peaks of cores MM2-9 shown in histogram, with Gaussian fittings overlaid 
in thick gray lines. 



Table 1. Measured parameters and masses of mm continuum cores 





R.A. 


Decl. 


Peak Flux 


Total Flux 


Mass a 




(J2000) 


(J2000) 


(Jybeam -1 ) 


(Jy) 


(M ) 


MM1 


23 h 13 m 45.358 s 


61°28'10.25" 


3.3 


3.6 


20.1 ± 11.6 


MM2 


23 h 13 m 44.724 s 


61°28'12.19" 


0.18 


0.37 


24.4 ± 12.7 


MM3 


23 h 13 m 46.432 s 


61 28'13.95" 


0.17 


0.31 


20.5 ± 10.7 


MM4 


23 h 13 m 43.729 s 


6i 28'11.35" 


0.11 


0.22 


14.5 ±7.8 


MM5 


23 h 13 m 45.496 s 


61°28'17.69" 


0.11 


0.30 


19.8 ± 10.6 


MM6 


23 h 13 m 46.942 s 


61°28'09.92" 


0.10 


0.15 


9.1 ±5.9 


MM7 


23 h 13 m 45.143 s 


61°27'57.95" 


0.082 


0.12 


5.2 ±2.9 


MM8 


23 h 13 m 43.446 s 


61°27'47.86" 


0.067 


0.10 


4.7 ±3.0 


MM9 


23 h 13 m 43.286 s 


61°27'59.49" 


0.049 


0.088 


6.2 ±3.0 



a Calculated with a dust opacity K225GHZ = 1 cm 2 g 1 , and temperatures of 245 
K for MM1, 40 K for MM2-5, and 43, 58, 54, 38 K for MM6-9; see the discussion 
in Section 13.11 for details about the mass uncertainties. 
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Table 2. Derived outflow parameters 



M out a 


P b 

-Tout 


^out 


^dyn 


Mout 


-Pout 


(M ) 


(M kms" 1 ) 


(erg) 


(y r ) 


(M yr- 1 ) 


(M Q kms' 1 yr" 1 ) 


50 


4.6 x 10 2 


4.9 x 10 46 


2 x 10 4 


2.5 x 1(T 3 


2.3 x 10-2 



a Calculated for outflowing gas at < —64 kms 1 and > —48 kms 1 ; 

b Without correcting for an unknown inclination angle; this affects all the pa- 
rameters except M out ; 
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Table 3. Molecular lines detected in MMl a 



rrequency 


Line c 


-E'up 


Peak d 
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(K) 


(Jy beam -1 ) 
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1 00 
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Zol 


l.OO 


ziy. /yoozu 


nwpn 1 n n 
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KQ 
DO 


1 £Q 
l.OO 
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c;n 

oy 


U.00 


oon 1 7Q 1 nf; 
ZZU.l ( oiyo 


pn 00 11 in 
Uri 2 bU lii,ii-iU ,io 


1 1 


n ns 


01 n nsQfi'7^ 


Uii 3 Uii Z0 3 ,22^4 4i 2o-^ 


snQ 


n ik 
u. ( 


01 n nnQfi^s 

ziy.yyoooo 


r^jj piu oq 00 z? 
UM3UII Z0519— ZZg,i7i!/ 


777 
III 


U.04 


oon ri7s^fii 
ZZU.U /oDOl 


(~<tt p>tr s it? 
Uri 3 Uii 00,8-' 1,6^ 


n7 

y < 


no 
z.yz 


oon fifif;7S/i 
ZZU.OOO ( 84 


pn fin 01 01 


1 1 so 
lloZ 


n /i k 


oon ^snn^fi 

zzy.ooyuoo 


pn pitr 1 c; 1 f; z? 
Uri 3 Uri 15 4 ,ii-lb 3 ,i 3 ^ 




/i c;n 
4.i>y 


oon vrcov^e 
ZZy. 100 1 00 
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c;7n 


1 nQ 
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OQn no7n/i7c 
Zoll.UZ / U4 / 


pn Q /I Z? 


/in 
4U 


Q Qfi 

0.0O 


230.368763 


f'Tr ATT OO Ol Z7 1 

OH3UII ZZ4,i8— ZI517.C/ 


684 


1 

1.56 


220.362774^ 


13 CH 3 OH 17_3.i5-i8o.i8 


410 


0.41 


220.475807 


CH 3 CN 12 8 -11 8 


526 


0.39 


220.539324 


CH 3 CN I27-II7 


419 


0.71 


220.594423 


CH 3 CN 12 6 -11 6 


326 


1.87 


220.641084 


CH 3 CN 12 5 -11 5 


248 


1.66 


220.679287 


CH 3 CN 12 4 -11 4 


183 


2.07 


220.709017 


CH 3 CN 12 3 -11 3 


133 


2.98 


220.730261 


CH 3 CN 12 2 -11 2 


98 


2.79 


220.743011 


CH 3 CN 12i-lli 


76 


3.05 
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Table 3 — Continued 



Frequency 13 


Line c 




Peak d 


(GHz) 




(K) 


(Jy beam -1 ) 


OOn 7/1 70£?1 

221). /4 / Zol 


r^u i o 11 


oy 


O 7C 

2. 10 


01 n 0^071 £ 
Ziy.Zoo / 10 


r^TJ r^TJC\ Q1 on 
Uil 3 UilU Ol5 27 - <J^2,30 


01/ 


n a q 
U.4o 


OQn Qm ssn 
ZoU.oUlooU 


Uri3UrlU 1^2,11-112,10 


fi 1 

ol 


n £/i 
U.04 


9 1 O /111 7fi Q 

Ziy.411 /Uo 


pit npun 1 fi 17 
l^li 3 UUliU 15io,8-1 'io,7 


OOO 


n qo 
u.oy 


Ziy. /0oy4 1 


r^tr or^tro is 1 7 
L/H3UL/HU 109,9—1^9,8 


Q A O 
o4Z 


n fi/i 
U.04 


01 o sooi of; 
Ziy.oZZlZO 


L/li 3 UOriU I5i 0l 9(8) _i ' 10,8(7) 


OOO 


n a 1 
U.41 


oon riQnQQO 
ZZU.UoUooy 


r^u npun 1 s 17 

OllsUollU log, 10— 1*9,9 


Q /I O 
o4Z 


n c;7 
u.o / 


oon 1 k^sss 

ZZU. 100000 


pit npnn 17 1 fi 

Ori3^J01HJ I1413— IO412 


1 no. 

1UO 


1 Ofi 
l.ZO 


oon i nnoc k 
ZZU.iyUZoO 


r^u npnn 17 it? 
Uil3UUllU 1/4,13—104,12 


1 no 
lUo 


1 QO 

l.oZ 


oon o^snQR 
zzu.zoouyo 


OH3UOHU lO8,10 _i '8,9 


001 


n 41 

U.41 


oon Qfins77 
ZZU.obyo / / 


Pitt npun 1 s 17 
Url 3 UbilU 1» 8 ,10^1'8,9 


ool 


n c;o 
U.oZ 


oon m QnKK 
zzU.y loyOD 


p'u npun 1 s 17 

OllsUollU 107,12— 1 1 7,11 


Q01 

oZl 


n Qf; 
U.00 


oon n/i£Qco 
ZZU.y4ooOZ 


OfUO 1 C 17 

oil3UoilU 107,11—1/7,10 


OOI 

621 


n o /i 
U.o4 


oon n77ns/i 
zzu.y / /yo4 
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U.00 
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ZZU.yoOooU 


p'tr nrun 1 s 17 
0H3U0IIU 107,11— 1 1 7,10 


qoi 
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n /i ^ 
U.40 


oon nnsQQc; 
ZZU.yyoooO 


UH3U0HU I8i5, 4 -1 i 15,3 


c;n 
Z0U 


n c;o 
U.oZ 


ooi n/i77m 
ZZ1.U4/ /yi 
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Ull3UbilU I«i4,5(4)-1 i 14,4(3) 
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PIT nPTTfl 1 fi 17 
UII3UUIIU IO145— 1 ' 14,4 


oqi 
261 


n Qfi 
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L-li 3 UL<ilU l«i3, 6 (5)-l < 13,5(4) 


Ol Q 
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n n7 

u.y / 


oon /inc;noi 
ZZy.4U0UZl 


P'u npun 1 fi 17 
L/H3U0IIU 103,15-1/3,14 


111 
111 


1 /I /I 
1.44 


ooq /ionQ/io 

ZZy .4ZU04Z 


pit npun is 17 
urisUL-nu 103,15-1/3,14 


111 
111 


1.40 


230.293951 


CH3OCHO 22g 13— 228,14 


204 


0.75 


230.315800 h 


CH3OCHO 22g, 14— 228,15 


204 


0.69 


230.878810 


CH3OCHO 18^14-1714,13 


301 


0.45 


220.601927 


C 2 H 5 OH 13i,i 3 -12 ,i2 


75 


1.08 


229.491131* 


C 2 H 5 OH 17 5 ,i2-174,i3 


160 


0.66 


229.818039 


c 2 H 5 OH 265,22-264,23 


385 


0.43 


230.177998 


C 2 H 5 OH 45 2 ,44-44 3 ,43 


894 


0.41 


230.230743 


C2H5OH 132,n-12 2 ,io 


143 


0.66 


230.473005 j 


C 2 H 5 OH 18 3 ,i 6 -174,i4 


215 


0.38 
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Table 3 — Continued 



Frequency 13 


Line 


E up 


Peak d 


{(jtlZ J 




(K) 


(Jy beam -1 ) 


230.672554 


C2H5OH 132,ii-12 2 ,io 


139 


0.66 


230.793764 k 


C 2 H 5 OH 65,1-54,1 


105 


0.53 


230.953778 


C 2 H 5 OH 16 5)11 -16 4) i2 


146 


0.75 


230.991374 


C 2 H 5 OH 14 0)14 -13 1)13 


86 


1.08 


220.893110 


CH3OCH3 234,2o-23 3 , 2 i 


275 


0.40 


230.141374 1 


CH3OCH3 254, 22 — 253 22 


319 


0.51 


229.234 m 






0.29 


230.578 m 






0.65 



a We list lines (excluding CO and 13 CO 2-1) detected above 
6a (0.21 Jybeam" 1 ); 

b Rest frequencies obtained from the JPL catalog at 
http : / /spec . j pi . nasa. gov /ftp/ pub / catalog / catf orm . html [ 

c Quantum numbers in parentheses denote degeneracy, for 
instance, HNCO 10 3)8 (7)-9 3) 7(6) represents HNCO 10 3 , 8 -9 3 ,7 
and HNCO 10 3 , 7 -9 3 ,6; 

d Derived with the MIRIAD task MAXFIT; 

e Under LTE, the peak flux is unexpectedly higher than 
other CH3OH lines (e.g., CH 3 OH (8_i, 8 -7 , 7 )£), which may 
be partly due to blending with C3H7CN (52 6)46 -5l6,45), 
^ O =230.027116 GHz and C 2 H 3 CN (4 3>1 -4 2i2 ), z/ o =230.026216 
GHz; 

f Blended with CH3OCH3 (28 5 , 2 3-2 7 6 , 22 ), z/ o =230.368178 
GHz; 

g Not listed in the JPL catalog, but included in Cologne 
Database for Molecular Spectroscopy; 
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h Blended with CH 3 CHO (12 2i i 1 -ll 2 ,io), z/ o =230.315740 
GHz; 

'Blended with CH 3 OCHO (31 5 , 27 -31 4 , 2 8), ^=229.492067 
GHz; 

JBlended with CH3OCH3 (10 8>3 -11 7 ,4), ^=230.473591 GHz; 

k Blended with C 2 H 5 OH (6 5)2 -5 4 , 2 ), v Q =230. 793864 GHz; 

'Probably blended with C 3 H 7 CN (36 6) 3i-35 5i3 o), 
z/ o =230. 142063 GHz; 

m Not identified. 



